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C^l , Relationships are obtained expressing the breaking of spin-reversal symmetry by an external 

magnetic field in Gibbsian canonical equilibrium states of spin systems under specific assumptions. 
These relationships include an exact fluctuation relation for the probability distribution of the 
magnetization, as well as a relation between the standard thermodynamic entropy, an associated 
spin-reversed entropy or coentropy, and the product of the average magnetization with the external 
field, as a non-negative Kullback-Leibler divergence. These symmetry relations are applied to the 
model of noninteracting spins, the ID and 2D Ising models, and the Curie- Weiss model, all in an 
external magnetic field. The results are drawn by analogy with similar relations obtained in the 
context of nonequilibrium physics. 
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I. INTRODUCTION 
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The breaking of a discrete symmetry by an external field is a common phenomenon in condensed matter at equi- 
librium. Examples are given by spin systems in a uniform external magnetic field, such as the Ising or Curie- Weiss 
models. In the presence of the external field, the Hamiltonian of these systems is mapped onto the Hamiltonian in the 
opposite field under spin reversal. Concomitantly, the system acquires a magnetization, which breaks the spin-reversal 

■ symmetry of the zero-field Hamiltonian. In finite systems or in finite regions of infinite systems at equilibrium, the 
order parameter given by the magnetization fluctuates around its average value according to some probability distri- 
bution and we may wonder if this probability distribution obeys general relationships resulting from the underlying 

i ^ i ' spin-reversal symmetry. 

The similarity is striking with the situation in nonequilibrium open systems in contact with reservoirs where the 

t-H ■ thermodynamic forces induce currents breaking the time-reversal symmetry. Nevertheless, the current fluctuations 
obey remarkable symmetry relations known under the name of fluctuation theorems [JQ. The question arises 
whether similar results can be obtained for the fluctuations of the order parameter in equilibrium systems where 
discrete symmetries such as the spin-reversal symmetry are broken. 

■ The purpose of the present paper is to show that it is indeed the case for equilibrium systems under specific 
• ' assumptions. We consider spin systems described by Gibbsian equilibrium canonical distributions. The Hamiltonian 

is supposed to depend linearly on the order parameter multiplied by the external field, as it is the case in the Ising and 
Curie- Weiss models. The fluctuations of the order parameter are characterized in terms of the equilibrium probability 
distribution of the magnetization and associated large-deviation functions such as the cumulant generating function 
of the magnetization. Using the aforementioned symmetry of the Hamiltonian, a fluctuation relation is here proved 
in a form that is similar to those previously obtained for nonequilibrium steady states. 

Furthermore, a spin-reversed entropy or coentropy is introduced besides the standard thermodynamic entropy in 
order to characterize the probability of spin configurations that are opposite to the typical configurations. The spin- 
reversed entropy is analogous to the time-reversed entropy that has been previously introduced in the context of 
nonequilibrium statistical mechanics |9j. Here, the difference between the spin- reversed and the standard entropies 
turns out to define a Kullback-Leibler divergence that is always non-negative and proportional to the external field 
multiplied by the average magnetization. This result completes the analogy with the time-reversal symmetry breaking 
by nonequilibrium steady states, for which similar relationships have been established [il-iioj. The theory is illustrated 
with several systems: (1) a system of noninteracting spins in an external magnetic field; (2) the one-dimensional Ising 
model; (3) the Curie- Weiss model; (4) the two-dimensional Ising model. The present paper extensively develops 
preliminary results reported in Ref . . 

The paper is organized as follows. The framework is presented in Section [Til The equilibrium fluctuation relation is 
proved in Section [Ml In Section ITVl the spin-reversed entropy is introduced and shown to combine with the standard 
entropy to get the non-negativity of the average magnetization multiplied by the external magnetic field. The theory 
is applied to the system of noninteracting spins in an external magnetic field in Section [Vl to the one-dimensional 



2 



Ising model in Section IVTI to the Curie- Weiss model in Section IVlIl and to the two-dimensional Ising model in Section 
IVIIII Conclusions are drawn in Section IIXI 

II. SPIN- REVERSAL SYMMETRY AND ITS BREAKING BY AN EXTERNAL MAGNETIC FIELD 

A. The Hamiltonian and its symmetries 

We consider systems composed of N spins a = {(Ji\f =l - For spin one-half, the individual spin variables take the 
values <7i = ±1 so that the state space S = {a} contains 2 N spin configurations. For spin one, the values would be 
di = +1,0, —1 and the state space would contain 3 N spin configurations. 

The energy of the system is given by the Hamiltonian function i/jv(<x; B) where B denotes the external magnetic 
field. We introduce the magnetization 

N 

M N (a) = Y,<ri (1) 

i=l 

which plays the role of order parameter. The Hamiltonian is assumed to depend linearly on the external magnetic 
field multiplied by the magnetization 

H N (a; B) = H N (a; 0) - B M N {a) (2) 

as it is the case for instance in the Ising or Curie- Weiss models. 
We consider the discrete symmetry of spin reversal: 

a R = Ra = -a (3) 

which is an involution because R 2 = 1. Accordingly, the spin-reversal transformation generates the discrete group 
Z2 = {1,.R}. We notice that the set of all the spin configurations {cr} that defines the state space is mapped onto 
itself by spin reversal, RT, = R{<r} — {—cr} = {cr} = £, because the state space contains the reversal of every spin 
configuration. 

In the absence of external magnetic field, the Hamiltonian is supposed to be symmetric under spin reversal while 
the magnetization is reversed: 

RH N (a;0)R = H N (a;0) (4) 
RM N {a)R = -M N {a) (5) 

As a consequence, spin reversal maps the Hamiltonian function in the external magnetic field B onto the Hamiltonian 
in the reversed magnetic field — B: 

RH N (a; B)R = H N (a; —B) (6) 

Therefore, the external field is expected to induce a magnetization, which breaks the Z2 symmetry of the Hamiltonian 
in the absence of external field. 

Moreover, the Hamiltonian i?jv(f;0) may also be symmetric under the transformations of another group G of 
transformations. This group depends on the geometry of the network of interactions between the spins. If the spins 
interact between the nearest neighbors of a lattice on a torus, the group G contains all the translations of the lattice, 
as in Ising models. If all the spins interact together in a fully connected graph, the group G is composed of all the N\ 
permutations of the vertices of the graph: G = SymTV, as in the Curie- Weiss model. 

B. Gibbs' canonical equilibrium states and thermodynamics 

The system is supposed to be in equilibrium at the temperature T, in which case the invariant probability distri- 
bution is given by Gibbs' canonical equilibrium state: 



(7) 
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where ft — (kT) 1 is the inverse temperature and k Boltzmann's constant. The normalization of the probability 
distribution ^2 ji^{a; B) = 1 is guaranteed by the partition function 

Z N (B) = tr e-f" 1 " (<T:S) = e ~ " N ( g ) 

Because of the symmetry (JHJ), the partition function is an even function of the magnetic field: 

Z N {B) = Z N {-B) (9) 
since summing over all the spin configurations is equivalent to summing over all their reversals: XLr(') = S-o-(') = 

£«,(■)■ 

The Helmholtz free energy is defined by 

F N (B) = -kT In Z N [B) (10) 
The energy is given by the statistical average of the Hamiltonian function 

E N = (H N ) B =trii N H N (11) 
with respect to the Gibbsian probability distribution ([7]), while the entropy is defined by its thermodynamic relation: 

S» = (12) 

so that the basic relation defining the Helmholtz free energy is satisfied: Fjf — En — TSm- 

The statistical average of the total magnetization is obtained by taking the derivative of the free energy with respect 
to the external field: 

(Mn)b = (13) 



The free energy per spin is defined by 



and the average magnetization per spin by 



f{B)= lim ±F N (B) (14) 

N— >oo ly 



Hb = ^(m n ) b = -^ (15) 

Now, the Legendre transform of the free energy may be introduced as 

g(m) = f(B)+mB (16) 

where B = B(m) is the magnetic field corresponding to the given magnetization m = — dB.f[B(m)]. Reciprocally, the 
free energy per spin can be recovered by the Legendre transform of the new function g(m), in which the magnetization 
per spin is replaced by its value corresponding to the external field B = d m g[m(B)\. In systems with phase transitions 
where these derivatives may not exist, the Legendre transform should be replaced by its extension called the Legendre- 
Fenchel transform [l^, [H[ ■ 

We notice that the symmetry ^ of the partition function results into the symmetry relations: 

f(B) = f(-B) and g{m) = g(-m) (17) 

Furthermore, the magnetic susceptibility per spin is defined as 

1 d(M N ) B 

XB = hm — — -- — (18) 

iV-yoo N OB 

and the variance of the magnetization fluctuations as 

a\ = lim i {(M 2 N ) B - (M N )%) (19) 
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III. SYMMETRY RELATIONS 
A. The fluctuation relation 

We introduce the probability Pb(M) that the fluctuating magnetization would take the value M = M^(cr) as 

P B (M) = {5 MMN[a) ) B (20) 

where 5m, w denotes the Kronecker delta function defined for M, M' £ Z and such that 5m, M' = 1 if M = M' and 
zero otherwise. This probability distribution is normalized to unity according to Pb(M) = 1. The statistical 
average is carried out over the Gibbsian canonical equilibrium distribution ([7]). 
The fluctuation relation is established as follows: 

Pb(M) = *—'£ e -l>*»W+f>B*»V , (21) 

- 1 V ~f3H N (a R ;Q) + f3BM N (a T< ) r zoo** 

= z^B) E e-^^-^W *m,-m kW (23) 

= E e-^("°)+^(*) 5_ M ,M N{a) (24) 

= e 20BM P B (-M) (25) 

From the first to the second line, the sum over all the spin configurations is equivalent to the sum over the spin- 
reversed configurations since they belong to the same set: i?S = S. From the second to the third line, the symmetries 
flU and (© are used so that H N (a R ;0) — H N (a;0) and M N (a R ) = — Mjv(<r). From the third to the fourth line, the 
Kronecker delta function allows us to restore the canonical equilibrium distribution by factorizing exp(2/3 BM) out of 
the sum. Finally, we obtain the fluctuation relation: 

Pb(M) _ 2pBM 



Pb(-M) 



C 2PBM (26) 



as an exact relationship under the aforementioned assumptions. This relation is one of the main results of this paper. 
We notice the analogy with the well-known fluctuation theorems for nonequilibrium steady states [1-8]. This analogy 
shows that these symmetry relations express the breaking of the discrete symmetry by the probability distribution in 
the presence of external constraints, here given by the magnetic field B. Indeed, in the absence of the magnetic field 
-6 = 0, the opposite fluctuations of the magnetization have equal probabilities and the symmetry is thus recovered. 
If the external field is non vanishing, a bias appears between opposite fluctuations and the fluctuations become more 
probable in one particular direction, which nicely expresses the breaking of the symmetry by the external field. 

B. The large-deviation function 

In the large-system limit, the probability distribution of the total magnetization is expected to behave exponentially 
in the number N of spins. In this regard, we can introduce the large-deviation function 

$n(m)= lim -— In P B {Nm) (27) 

TV— >oc is 



Therefore, the probability distribution (f20[) should behave as 

P B (M)~e- OTB W*! for N^oo (28) 



up to a subexponential dependence on N [12l - ll4j . 

Now, the fluctuation relation ([2"o]) implies that the large-deviation function (j2~?| obeys the following symmetry 
relation: 



<& B (-rn) - <&b(to) = 2j3Bm 



(29) 
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as can be verified by a straightforward calculation. 

The large-deviation function can be expressed in terms of the Hclmholtz free energy per spin and its Legendre 
transform (fTB]) according to 

$ B (m)=p\g(m)-Bm-f(B)] (30) 
In order to obtain this result, we notice that the partition function can be written as 

Z N (B) = £ e ^M£ e -^(„ ; o) (31) 

M a 

by introducing the Kronecker delta function. Since the partition function behaves as Zn(B) ~ e~ N Pf( B ) in the 
large-system limit, we see that the sum over the spin configurations should behave as 

J2 e -f>»H*;o) §NmtMNi<r)r ^ e -N Mm ) for (32) 

in order to recover the function (fTB]) as the Legendre transform of the free energy per spin f(B). Since the probability 
distribution of the magnetization can be obtained from the expression 

Pb{M) = -\- eP»" J2 e-'*"W> 6 M;MN{a) (33) 

we find that the large-deviation function is indeed given by Eq. (|30[) in the large-system limit. 

We notice that the symmetry relation (f29| is verified because of the symmetry (|17|) of the function g(m). 



C. The cumulant generating function and its symmetry 

In the framework of large-deviation theory [l3l . [l4| , the cumulant generating function of the magnetization fluctu- 
ations is defined as 

Q B (X) = lim -1 ln(e- AM ~) B (34) 

N—>oo iv 

This function is generating in the sense that the average magnetization per spin, its variance, as well as all the 
higher-order cumulants can be obtained by taking successive derivatives with respect to the parameter A: 

, , 9Q B , . 2 9 2 Q B , , , . 

( m )B = -gT"(°) » °"B = dj^(°> ' - ( 35 ) 

Now, the statistical average in Eq. (I3"4l can be written in terms of the probability distribution ([20)1 of the magne- 
tization as 



. B = ^P B (M)e- AM (36) 



M 



Because of the fluctuation relation (|26p. we immediately obtain the symmetry relation of the cumulant generating 
function: 

QbW = Q B (2j3B-X) (37) 

which is similar to symmetry relations obtained elsewhere for nonequilibrium systems 

In analogy with the results of Ref. [ 1 51 ] , the response coefficients can be related to the cumulants of the magnetization 
fluctuations be taking successive derivatives of the symmetry relation (|37l) with respect to the parameter A and the 
external field B. In particular, taking one derivative d\ and another &b gives 

and setting A = B = 0, we recover the well-known relation 

Xo=P<t 2 (39) 
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between the magnetic susceptibility and the variance of the magnetization fluctuations in the absence of external field. 
Similar relationships can be obtained between higher-order response coefficients and cumulants. We notice that, for 
Gibbsian equilibrium states, the relation (|39| also holds for non- vanishing external fields, xb = Pa%, which is not 
the case in nonequilibrium steady states since these latter do not have the simple Gibbsian form ([7|). 

The cumulant generating function can also be obtained in terms of the partition function and the free energy per 
spin. Indeed, we have that 

(e- AM "> B = Y m {a- B) e - XM -^ = V e -/UM»iO)+(/«-A)JM») = Zjv( f ~f" lA) (40) 

^ Z N (B) ^ Z N (B) 

In the thermodynamic limit N — ¥ oo, we thus obtain 

Q B (X)=/3[f(B-fi- 1 X)- f(B)] (41) 

In this regard, the symmetry for the partition function or (|17[) for the free energy per spin implies the symmetry 
relation (|37[) for the cumulant generating function. 

A further consequence of Eq. ([36]) is that the cumulant generating function ([34]) can be obtained as the Legendre 
transform of the large-deviation function (|27l) : 

Qb(X) = $s[m(A)] + Am(A) where [ m (A)] = -A (42) 

am 

Conversely, the large-deviation function is given by the Legendre transform of the cumulant generating function: 

qq b 

&B(m) — Qs[A(m)] — rn\(m) where [A(ra)] = m (43) 

OA 



The Legendre-Fenchel transform should be used if these functions were non differentiable [12Ml 



IV. ENTROPY, COENTROPY AND BROKEN SYMMETRY 



In this section, we complete the parallelism with the results obtained for nonequilibrium steady states by considering 
a spin-reversed entropy or coentropy defined in analogy with the time-reversed entropy per unit time that has been 
previously introduced [9]. This quantity is shown to combine with the standard thermodynamic entropy to form a 
Kullback-Leibler divergence, as well as an important relationship with the average value of the magnetization. 



A. Entropy and coentropy 

The standard thermodynamic entropy is defined by Eq. (|12p . from which we get the usual expression: 

Sn = — fctr/ZArln/xjv (44) 

for the Gibbsian equilibrium state J?]). 

The breaking of the spin-reversal symmetry by the external field can be characterized by comparing the probability 
/iAr(cr; B) of some spin configuration a with the probability of the reversed configuration: 

l4f{a;B) = fi N (Ra;B)=n N (-a;B) (45) 

Here, we introduce the spin-reversed entropy 

S% = -fctr MJV ln^ (46) 

in analogy with a similar quantity introduced elsewhere in the context of nonequilibrium statistical mechanics Q . 
Now, the difference between the spin-reversed and the standard entropy defines the Kullback-Leibler divergence: 

D((i N \\^)=tifi N \n^. = l(S*-S N )>0 (47) 

which is known to be always non- negative [T6| . In this regard, the spin-reversed entropy could be called a coentropy 
since it combines with the entropy to form the non-negative Kullback-Leibler divergence. 
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The Kullback-Leibler divergence vanishes in the absence of external field when the probabilities of every spin 
configuration and its reversal are equal Hn = A* at- However, in the presence of an external field, the probability of 
a reversed configuration Rtr is expected to take a different value than the probability of the configuration a itself, 
in which case the Kullback-Leibler divergence becomes positive. Accordingly, the Kullback-Leibler divergence is a 
measure of the breaking of the symmetry at the level of the probability distribution /i^r. This is confirmed by the 
remarkable relation 

D(fM N \\f4) = I {Sn ~ S N ) = 2/35 (M N ) B > (48) 
which is the consequence of the fact that 

ixfi N hx^ = 28 B tv a N M N (49) 

holds for the Gibbsian equilibrium state (J7]) . A corollary of this result is that spin systems described by a Hamiltonian 
function ([2]) with the symmetries (j4j and ([5]) is always paramagnetic because the average value of the magnetization 
should point in the same direction as the external magnetic field. 
If we introduce the entropy and coentropy per spin as 

lim — S N (50) 



AT->-oo N 
1 



s R = Jim —S% (51) 



the relation (l48t writes 



I ( S R - s ) = 2^B (m) B > (52) 

where (m)s is the average magnetization per spin (|15|) . Consequently, we have the general inequalities: 

s R > s > (53) 
The relation (|52p is another of the main results of this paper. 



B. The disorders of the spin configurations and their reversals 

In order to interpret the previous result, we consider spin systems on an infinite c?-dimensional lattice Z d . In the 
infinite-system limit, the equilibrium states 

/i = lim /itjv (54) 

N^oo 

can be constructed as the so-called Dobrushin-Lanford-Ruelle states 0,11 Let A C Z d be a part of the lattice 
and |A| its volume. If a G is a spin configuration of the infinite lattice, (7\ is the spin configuration restricted on 
the domain A of the lattice. In this framework, the entropy and coentropy per spin are defined as 

k 

s = lim -y-r-r tr /z(<r A ) In /z(<r A ) (55) 

|A|->oo |A| 

s R = lim - — txfi(a A )hifi(a^) (56) 

|A|— >oo |A| 

These quantities characterize the disorders in the spin configurations and their reversals. Indeed, for almost every 
spin configuration a with respect to the probability measure fj,, the probabilities that the spins in the domain A are 
found in the configuration a a and their reversals in the configuration <j R should decay in the infinite- volume limit 
|A| — > oo as 

MaA^c-l^ (57) 
/.(O-e-W 5 "/' (58) 
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The faster the decay, the rarer the spin configuration a or its reversal a among all the possible configurations. 
Because of Eq. (1521 , we have the general relationship 

2 

s R = s + — B (m) B with B(m) B >0 (59) 

between the entropies (|55p - (|56[) and the average magnetization per spin (m)s. As a consequence, the ratio of the 
probabilities of opposite spin configurations behaves as 

MK) ^ e 2^(m> B |A| (60) 
A*(°a) 

for /i- almost every configuration a in the limit |A| — > oo. If the external field vanishes B — 0, the spin configurations 
and their reversals are thus equiprobable. However, in the presence of an external field breaking the symmetry, the spin 
reversal of every typical configuration is always less probable than the typical configuration itself because s R > s > 
for B 0. In this case, the spin-reversed configurations appear more disordered than the typical confi gur ations 
themselves. The analogy is here complete with similar results obtained for nonequilibrium steady states }9l llOj. 

We notice that there is an important difference between Eq. and the fluctuation relation (j2"6l . Indeed, Eq. (joT))) 
concerns typical spin configurations a characterized by the average magnetization per spin (m) b , although the fluc- 
tuation relation (|26p depends on the variable magnetization M and not on its average value (M)b- 

All these results are illustrated with specific models in the following sections. 



V. THE NONINTERACTING SPIN MODEL 



A. Hamiltonian and thermodynamics 



In this section, we consider the simple model of N noninteracting spins <Tj = ±1 in an external magnetic field B. 
The Hamiltonian of this system writes 

N 

H N (*;B) = -Bj2<n (61) 

i=l 

that is symmetric under the group G — SymiV and also under the group Z2 = {1, R} if B = 0. This latter symmetry 
is broken if B 7^ 0. 

For an equilibrium spin system at the temperature T = (fc/3) -1 , the partition function is given by 

Z N (B) = tre"^" = (2 cosh (3B) N (62) 

the free energy by 

F N = —NkT In (2 cosh (3B) (63) 

and the energy by 

E N = (H N ) B = -B(M N ) B (64) 

with the average total magnetization 

(M N ) B = N tanh/3B (65) 

The entropy of this system has Shannon's form 

Sn = Nk (— p + lnp + — p- lnp_) (66) 

in terms of the probabilities 

e ±/3B 

P± = c+ 0B +c -0B ( 67 ) 

that a spin would be up or down. These probabilities satisfy the normalization condition p + + p- = 1 . The average 
magnetization per spin is equal to (m)s = p+ — p- = tanh/3_B. 
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B. Fluctuation relation 

The probability (|20| that the fluctuating magnetization would take the value M can be exactly calculated for the 
noninteracting spin model. Indeed, this probability is given by the binomial distribution as 

ftw-s^i*;^- (») 

with the probabilities (|67[) and the numbers N± of up and down spins. These numbers are related to the value M of 
the magnetization according to 

N± = i (N ± M) (69) 



so that the probability distribution (|6"8]l can be written as 



jyl N+M N - M 



Pb{M) _ {P+\ M _ :,; W 



Therefore, we see that the fluctuation relation (1261) holds 

>. M 

P B (-M) \p_J [<1) 

because of Eqs. (J6TJ) . 

Here, the cumulant generating function (|34[) has the expression 

/ \ \ , cosh BB .„„, 

Qb (A) = In J" 72 
cosh(/3ii — A) 

so that the symmetry relation ([37"1) is satisfied. The Legendre transform (|43|) gives the large-deviation function 

$ B ( m ) = I±HL In + In - /35m + In (2 cosh BB) (73) 

This function determines the asymptotic behavior of the probability distribution (|70[) in the large-system limit N — > oo 
and its symmetry relation (|29[) is indeed satisfied. 

C. Statistics of typical and reversed microstates 

The probability of a microstate or spin configuration a with N± spins that are up or down is given by 

t i N {a) = P N + + P N - (74) 
and the probability of the corresponding spin-reversed configuration by 

M" R )=P+~i>- + (75) 

The ratio of both probabilities takes the exact value 

t \ / \ N+-N- 
J*W = (P±) =e 20B i N + -N.) (76) 

/M°" ) \P-J 

Now, let us consider these expressions for the most probable configurations, i.e., for the typical microstates such 
that 

N+~Np+ and N- ~ N P - (77) 

in the large-system limit N — > oo. Using the Stirling formula, the total number of these typical microstates can be 
evaluated as 

V ' - 1 „ e^ (78) 



JV+UVJ 1V ^°° ^Trp+p-N 
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in terms of the standard thermodynamic entropy per spin: 

s = —k (p + hip + + p_ lnp_) 



(79) 



If this entropy is non vanishing, the probability that one among all the typical microstates is observed is thus decaying 
exponentially as 



UN (a) 



-N- 



e N(p+ In p++p- In p_) _ e -Ns/k 



(80) 



at the rate given by the entropy per spin, while the probability of the corresponding spin-reversed microstates decays 
as 

-Ns R /k 



Hn(& R ) —N- 



at the rate 



e N(p- lnp + +p + lnp_) _ e 



-k (p- lnp + + p + hip J) 



(81) 



(82) 



that is the coentropy per spin. Here, we can directly verify that the ratio of both probabilities indeed satisfies the 
relation (|60p because, for /Zjv-almost all the microstates a, we have that 



1 AT( s R -s)/fc 



with 



^ (s R -s) = (p+ -p-) In — = 2/3Btanh^ J B = 2/3B{m) B > 



P- 



(83) 



(84) 



This quantity is non negative because it is related to the Kullback-Leibler divergence (|48p . Another way to verify the 
result (|60p is to use directly Eq. (|75|) on the subset of the most probable microstates for which we have that 



N+ — N- ~ 7V( p+ - p_) = AT tanh fiB = N (m) B 
leading to the same result ([55} with Eq. ([5i|l. 

all the microstates 



(85) 



typical microstates 



reversed microstates 



FIG. 1: Schematic representation of the set of all the 2 possible microstates containing the subset of the typical microstates 
with N+ = Np+ spins up and iV_ = Np- spins down together with the subset of the corresponding spin-reversed microstates 
with = Np- spins up and N+ = Np+ spins down. The grey area indicates the region where the probability distribution is 
concentrated. The probability distribution culminates on the typical microstates that are the most probable microstates. The 
subsets of typical and reversed microstates coincide if 7V+ = iV_ , but they do not overlap as soon as JV+ / N- . 



In the case of a non-vanishing external field B ^ breaking the Z2 symmetry, the situation is illustrated in 
Fig. [l] which depicts the set of all the possible microstates of the spin system. On this set, the Gibbsian equilibrium 
probability distribution is concentrated in a certain area shown in grey where we find the subset of the most probable 
microstates, called the typical microstates, for which Eqs. (fTTf hold. Because N + ^ iV_ if B ^ 0, the corresponding 
spin-reversed microstates form a disjoint subset in this case. The spin reversals of the typical microstates are thus 
found in an area of low probability because of the general relation (|83l) . 

In the model of noninteracting spins, the validity of the fluctuation relation and the other symmetry relations can 
thus be verified by direct calculations. 
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VI. THE ID ISING MODEL 
A. Hamiltonian and thermodynamics 

A further model where the symmetry relations can be verified analytically is the ID Ising model of Hamiltonian 

N N 

H N {a;B) = -J^a i a i+1 -B^(j i (86) 

i=l i=l 

where the N spins o\ = ±1 form a ring with <jn+i = fl- 
it is well known that the canonical partition function of this model can be calculated thanks to the transfer-matrix 
technique [20| . Accordingly, all the thermodynamic quantities are exactly calculable for this ID model. 



B. Fluctuation relation 



Here, we are interested in the generating function of the statistical moments of the fluctuating total magnetiza- 
tion M^((j) 



i r N N j_ 

(c- XMn )b = E ex P P J E WW + IP** ~ A ) E U - 1± Y ±1 

N ^ ' a |_ »=1 i=l 

which can be expressed as 

\rV N 

(.-»->- ^ 

in terms of the following transfer matrix 



(87) 



v * = ( e -£J e /JJ-,3s+A J (89) 

This transfer matrix is real and symmetric. Therefore, it has two real eigenvalues A + (£> — /3 _1 A) > A_(£> — (3~ 1 \). 
The trace of its ./V th iterate can thus be written as 

tr = A+(B - /3" 1 A) JV + A_(B - 0-^)* (90) 

Accordingly, the trace is dominated by the largest eigenvalue A+(B — f3~ 1 \) in the large-system limit N — > oo. 
Therefore, the cumulant generating function (|34|) is obtained as 



A+^-^A) 



»(A)=ln A ,„ +y ' (91) 



With the explicit expression of the leading eigenvalue A + , we find that 



cosh/?i?+ y/sinh 2 /3i3 + c-W 
Q B (A)=ln (92) 

cosh(/3B - A) + yjsinh 2 (/3B - A) + c- 4 ^ 

The symmetry relation (|37p is thus verified for this model as well. This symmetry is illustrated in Fig. [5] depicting 
the cumulant generating function for different values of (3J and f3B. In every cases, the function is symmetric under 
the transformation A — ► 2/3 B — A. This is obvious if B = 0. If (3B = 0.2 for instance, the symmetry holds under the 
map A — > 0.4 — A, as indeed observed. Another symmetry is also visible in Fig. [21 namely Qb(A) = Q-b(— A), which 
results from the symmetry of the partition function. 

Moreover, the symmetry A — > 2/3B — A already manifests itself at the level of the transfer matrix (|89| that satisfies 
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FIG. 2: Cumulant generating function (|92|) of the ID Ising model versus the generating parameter A for /3J = 0.5 (solid lines) 
and PJ = —0.5 (dashed lines). 




FIG. 3: The large-deviation function (|95[1 of the ID Ising model versus the magnetization variable m for the same values of 
the physical parameters as in Fig. [2] 



is the matrix performing spin reversal: <x; = ±1 — s- erf = =pl. Consequently, the symmetry A — > 2/3 B — A holds not 
only for the cumulant generating function determined by the leading eigenvalue A+, but also for both eigenvalues 
A ± . 

From Eq. (|92j) . the average magnetization per spin is given by 

/ \ d ® B m\ sinh/3B 

(m)B = -^r-(O) = , (94) 
Vsinh 2 JB + e~W 

as expected. The magnetization behaves smoothly in the plane (T, B) except at the point (T = 0, B = 0) where the 
function is discontinuous, although this ID model does not have a genuine phase transition at a positive temperature. 
The Legendre transform of the cumulant generating function (|9"!?|) takes the following form: 



VT -to 2 (cosh (3 B + Vsinh 2 (3B + e'^A 



\ yj 1 — ?7i 2 
$ b (W) = In m In — BB to (95) 

-20J + ^ m 2 e -4/3J + 1 _ TO 2 me -2/3J + ^ m 2 g-4/3,7 + J _ m 2 
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This function is depicted in Fig. [3] for the same conditions as in Fig. [2j The analytical expression (|95|) is such that 
the symmetry relation is indeed satisfied, hence the fluctuation relation (J2U) holds for the ID Ising model. 



C. Entropy and coentropy per spin 



The free energy per spin is given by 

/ = -kT \nA + (B) = -kT In ^cosh^B + \J sinh 2 (3B + e-^ J 
allowing us to obtain the standard entropy per spin by the usual thermodynamic formula 

dT 



(96) 



(97) 



Now, the coentropy per spin s R is calculated with Eq. (|59[) and the average magnetization per spin (|94[) . The entropy 
and the coentropy are depicted in Fig. 0] as a function of the external magnetic field B for different values of the 
temperature. Both entropies tend to increase with the temperature, as expected since they both characterize some 
form of disorder. If B — 0, the coentropy is equal to the entropy. If B ^ 0, the coentropy is larger than the entropy 
because of the symmetry breaking and the general inequality (|52|) . For \B\ — > oo, the entropy decreases because of 
the ordering of the spins in the direction of the average magnetization, while the coentropy increases because the 
spin-reversed configurations become rarer and rarer. The difference between the coentropy and the entropy behaves 
as 

s R - s ~ 2 k e 2fiJ p 2 B 2 for (3B < 1 (98) 
which explains that the separation between them is more pronounced at lower temperatures, as seen in Fig. @] 

1 



at 



0.8 
0.6 
0.4 
0.2 



k'l-5 \ \ 


i.i i r i i i i 


kT- 1 \ / 




- kT=0.5 ~" 

1 1 r"^ , i , ; , 1 



-0.5 




B 



0.5 



FIG. 4: The entropy (solid lines) and coentropy (dashed lines) per spin of the ID Ising model for J = 0.5 versus the external 
magnetic field B. k denotes Boltzmann's constant and T the temperature. We notice that the entropy per spin tends to 
s = k In 2 as kT > J. 



VII. THE CURIE- WEISS MODEL 



This model features a paramagnetic-ferromagnetic phase transition, which occurs at a positive critical temperature. 
Below the critical temperature, a phenomenon of spontaneous symmetry breaking induces a non- vanishing magnetiza- 
tion in the absence of external magnetic field. Therefore, the probability distribution of the magnetization is expected 
to be bimodal in finite systems, which is an interesting situation for the fluctuation relation. 
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A. Hamiltonian and thermodynamics 



The Hamiltonian function of the Curie- Weiss model can be written as 



H N (a;B) 



J 

2N 



M N (a) 2 - B M N {a) 



(99) 



in terms of the magnetization ([T]) (l2| . This Hamiltonian is invariant under the whole symmetric group G = SymiV 
because the Hamiltonian only depends on the total magnetization that has this invariance. The Hamiltonian is also 
invariant under the group Z2 = {1, R] if B = 0, but this symmetry is broken if B 7^ 0. 

The total magnetization takes the N + 1 different values M = N — 2n with n = 0, 1, 2, N. The number of spin 
configurations with the given magnetization M is equal to 



Cm 



AH 



' N+M \ 1 ( N-AI \\ 



2 )• \ 2 

In the large-system limit, this number can be approximated as 



C 



M 



irN{l - m 2 ) 



exp 



-N 



1 



m 1 
— In — 



1 — m 1 
— - — In — 



(100) 



M 

with to = — (101) 



by using Stirling's formula. 

Since the increment of the magnetization per spin is equal to dm = 2/N, the partition function can thus be written 

as 



+ 1 



dm ■ 



-Nip B (m) 



vT 



with the function 



ipB(m) = 



1 



TO 1 

In ■ 



1 



TO 1 

In ■ 



TO /3J 2 



m — /3Bm = Pg(m) — (3Bm 



2 2 2 2 2 

The integral can be performed with the method of steepest descents by expanding this function as 

ipB(m) = ipB(nii) + ^s(wj) (to - m^ + - i/)g(mi) (to - mi) 2 H 

= 



(102) 



(103) 



(104) 



around its stationary points {mi} such that ip' B (mi) — and ip'g(mi) > 0. The stationarity ip' B (mi) = selects the 
value of the magnetization per spin satisfying the self-consistent condition: 



mi = tanh(/3Jrni + j3B) 



(105) 



Above the critical temperature kT c = J, this equation has only one solution. Below the critical temperature, there 
are three solutions, two of which satisfy the condition of local stability ipg(mi) > 0. In the thermodynamic limit 
N — > 00, the partition function can thus be evaluated as 



Z N (B) ~ £ 



-Nip B (mi) 



y/1 - /?J(1 - mj) 



(106) 



where the sum extends over the locally stable solutions with ip'g(mi) > 0. Since the partition function is related to 
the free energy per spin by Zj^(B) ~ e~ N Pf( B ) ; this latter is given by 



(3f(B) = Min m {^ B (m)} = (3Mm m {g(m) - Bm} 



(107) 



Above the critical temperature, there is only one solution to Eq. (|105j) . but two should be considered in the 
expression (|106|) of the partition function below the critical temperature. Since both terms are decaying exponentially, 
the convergence should not be expected before the size is large enough: N 3> \tpB(mi) — ^3(^2) | _1 . 

Once the free energy per spin is calculated, the average magnetization per spin can be obtained by Eq. (|15[) . which is 
depicted in Fig. [5]for different values of the external field B. The critical temperature is here equal to kT c = J = 0.5. 
This figure allows us to appreciate the effect of the external field on the average magnetization per spin. In order to 
have a moderate symmetry breaking effect, the value B — 0.001 of the external field will be used in the following. 
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FIG. 5: The non-negative average total magnetization versus the thermal energy kT in the Curie- Weiss model for J = 0.5 and 
several values of the external magnetic field B. 



B. Fluctuation relation 



The probability distribution of the magnetization M is given in the Curie- Weiss model by 

Pb(M) = ^ C M ex P (!» M 2 + PBM) (108) 

where Cm is the corresponding number of spin configurations (|100p . This distribution has been computed for J = 0.5, 
B = 0.001, N — 100, and different values of the temperature across the phase transition, as depicted in Fig. [6] We 
observe that, below the critical temperature kT c = J = 0.5, the distribution is bimodal with two peaks centered around 
the two most probable values for the magnetization. Above the transition, the distribution becomes unimodal. The 
effect of the small external field B = 0.001 is to induce an asymmetry in the distribution toward the direction of the 
external field. Therefore, the average value of the magnetization has the sign of the external field. 

In order to verify the fluctuation relation (p2rj)l . the distribution Pb(M) is compared in Fig. |5]with the prediction 
that it should coincide with Pb(— M) exp(2/3BM). We see the nice agreement between both even below the critical 
temperature where the distribution is bimodal and quite different from a simple Gaussian distribution. 



C. Large-deviation function 

For a better understanding of the behavior of the probability distribution Pb(M) across the phase transition, we 
consider its approximation (1281) in terms of the large-deviation function (|27p in the thermodynamic limit N — > oo. 
Using the approximation (jlOll) . the expression (|108[) can be written as 



p ^ M ^^)i^mh^)^* ,M mth ■»-# < 109 > 

in terms of the function (|103[) . Comparing with Eq. (|28p , we recover the expression of the large-deviation function 
from Eq. (fTUgj) . 

The probability distribution (| 108[) for J = 0.5, B = 0.001, and N = 50, 100, 150 is compared with its large-deviation 
approximation (|109[) in Fig.[7]below and above the critical temperature. In this approximation, the partition function 
is evaluated with Eq. (I106p . Above the critical temperature, the distribution is unimodal and a single peak contributes 
to the partition function. In contrast, two peaks contribute to the partition function below the critical temperature 
at kT = 0.4 in Fig. [7^. Still for N = 150, the two terms of Eq. (jl06[) are required to get the agreement seen in Fig. [7^,. 
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FIG. 6: The probability distribution Pg(M) of the magnetization M (open squares and lines) versus the magnetization M in 
the Curie- Weiss model for J = 0.5, B = 0.001, N = 100, and: (a) T = 0.3; (b) T = 0.4; (c) T = 0.5; (d) T = 0.6. The crosses 
show the values of Pb{— M) exp(2/3BM). The coincidence of the crosses with the squares is the prediction of the fluctuation 
relation J26J). 



Indeed, the dominant peak only contributes to 63 % of the partition function for N = 150. The reason is that the free 
energy per spin is given by /3f(B) = ipB{rni) = —0.7308, while the value corresponding to the subdominant peak is 
equal to ^b(to2) = —0.7272. Therefore, the second peak only becomes negligible if JV |^s(mi) — 4>B{m2)\~ 1 — 278 
for significantly larger sizes than TV = 150. 



D. The cumulant generating function 



The cumulant generating function Q34[) can be obtained using the large-deviation function f|103[) . First, the gener- 
ating funciton of the statistical moments can be evaluated as 

1 „ _ p -JV[*jj(m,)+Am,] 

<e- AM ~) B - £ h Wm (U0) 

Z N {B) ^ - - m f) 



17 



0.1 



0.01 



(a) 



0.001 - 



0.0001 







: ■ 


I Vl \ /V = 50 / ft \ 

ft \\ ^iliWi^ // \ 




■ I 


r \ \ 7v = ioo / / 1 






\n = 15o/ 












-1.2 -0.8 -0.4 0.4 
m = M/N 



0.S 



1.2 



0.1 



0.01 



0.001 r 



0.0001 



(b) 



TV = 50 .^ ssw 



. N = 100 



-1.2 




-0.4 0.4 
m = M/N 



1.2 



FIG. 7: The probability distribution Pb(M) of the magnetization M versus the magnetization per spin m = M/N in the 
Curie- Weiss model for J = 0.5, B = 0.001, and N = 50, 100, 150 (filled squares, circles, and diamonds) compared with 
its expression (|109[) in terms of the large-deviation function (|103|) in the large-system limit (lines): (a) below the critical 
temperature at T = 0.4 where the probability distribution is bimodal; (b) above the critical temperature at T = 0.6 where the 
probability distribution is unimodal. 
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FIG. 8: The cumulant generating function Qs(A) divided by the average magnetization per spin (m)s versus the rescaled 
generating parameter \/(2j3B) in the Curie- Weiss model for J = 0.5, B — 0.001, and different values of the temperature T 
across the phase transition at kT c — J — 0.5. 



where the sum extends over the roots rhi of ^>' B {fhi) + A = with ^(m;) > 0. These roots are here the solutions of 

rhi = tanh(/3 Jmj + f3B — A) (111) 
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Since the partition function behaves as Zn(B) ~ e Ar ^^( s ) and (e XMn ) b ~ e N Q B W in the large-system limit, the 
cumulant generating function is obtained as 

Q B (X) = Min m {i/> B (m) + Am - 0/(5)} = Min m {$ B (m) + Am} (112) 

which is valid either below or above the phase transition. 

For J = 0.5 and B = 0.001, the cumulant generating function has been calculated by Eq. (|112p and is shown in 
Fig. [8] for different values of the temperature across the phase transition after rescaling by the average value of the 
magnetization per spin. We observe that, above the critical temperature kT c = J = 0.5, the cumulant generating 
function is differentiable where the magnetization distribution is unimodal. In contrast, the function becomes non- 
diffcrcntiable below the transition where the magnetization distribution is bimodal. The reason is the presence of two 
competing solutions of Eq. (jllip as the parameter A varies. The minimum is exchanged between these two solutions at 
the critical value A = f3B, which creates the tent-like shape of the cumulant generating function. It should be noticed 
that the cumulant generating function is never exactly piecewise linear in spite of the appearance. Its nonlinearity 
tends to increase with the value of the external field B. 

The convergence is very slow near the non-differentiable point in the cumulant generating function. By a similar 
reasoning as for the partition function, we should expect convergence for N ^> \if) ^(mi) + Ami — V'sC 7 ^) — A7712I 
where rhi — rhi(X) are the two roots of Eq. (Illip . However, this lower bound on N diverges as |A — /3-B| -1 near the 
point of non-differentiability where the convergence is thus arbitrarily slow. 

In any case, the symmetry relation (|37[) is always satisfied above or below the phase transition. We notice that, in 
the ferromagnetic phase, both peaks of the magnetization distribution must be taken into account in order to obtain 
the symmetry (|37[) of the generating function. The symmetry would not hold if only one peak would be considered. 

VIII. THE 2D ISING MODEL 

A. Hamiltonian and canonical equilibrium states 

The last example studied in this paper is the famous 2D Ising model on a square lattice with nearest-neighbor 
interactions [ill [2(|. For finite systems, periodic boundary conditions are taken so that the square lattice forms a 
torus of size L in both directions, containing N = L 2 spins one-half. The Hamiltonian of this system is given by 

L L 

H N (a;B) = [-J cn,j +<r%,i+r) - Ba iyj ] (113) 

i=i j=i 

with CTi.L+i = o"i,i and OL+i.j = cij for i,j = 1,2,...,L. This Hamiltonian is invariant under the group G of 
symmetries of the toral square lattice. The further Z2 symmetry holds if B = 0, but is broken if B ^ 0, as in the 
previous examples. 

The 2D Ising model is exactly solved in the absence of external field for B = (l9l . [2pj . Otherwise, Monte Carlo 
tec hniq ues are available for computations in the presence of an external field, such as the algorithm by Metropolis et al. 
[2ll |22| . we here use. 

The 2D Ising model satisfies the symmetry conditions ((H) and (JS|) required for the fluctuation relation (121)1) , as well 
as for the relation (|4"5|) between the entropy and the coentropy. 

B. Fluctuation relation 

Here, we consider in detail the fluctuation relation (|26[) for the probability distribution of the magnetization com- 
puted by Monte Carlo simulations at different temperatures, J = 1, and in the presence of an external magnetic field 
B = 0.01 breaking the Z2 symmetry. We compare two lattices of finite sizes L = 10 and L = 20. We notice that the 
paramagnetic-ferromagnetic pha se transition at the critical temperature kT c = 2.269 J is an exact property only in the 
infinite-system limit L — > 00 |19j |. For systems of finite size, the transition manifests itself only as a crossover. At low 
temperature, the magnetization distribution is bimodal if the external field is not too large and unimodal otherwise. 

This distribution is obtained by sampling the magnetization every 10 x L 2 Monte-Carlo spin flips generated with the 
algorithm by Metropolis et al. [2l|,[22|]. The histogram of the magnetization is obtained with a sample of 10 7 values 
for the lattice of size L — 10 and 5 x 10 7 values for the size L = 20. The histogram gives the probability density 

PB(m) = -±- Pb{ - M ) ( 114 ) 

Am ^— ' 

N(m-Am/2)<M<N(m+Am/2) 
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where m = kAm with k = 0, ±1, ±2, ±K, Am = 1/K, and K = L. 

(a) 



(b) 




10 



sq 

0.1 



0.01 



L = 20 



-1.2 -0.8 -0.4 




m 



0.4 0.8 1.2 



FIG. 9: The magnetization histogram (squares and lines) for the 2D Ising model with J = 1, B = 0.01, and kT = 2.2 on a 
square lattice of size: (a) L = 10; (b) L = 20. The crosses show the prediction of the fluctuation relation (|115[) . 



In Fig. \§\ the histogram (|114l) of the magnetization is shown for lattices of sizes L — 10 and L = 20 below the 
transition at kT c = 2.269. As the size increases, we see that the bimodality of the magnetization distribution becomes 
stronger. 
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FIG. 10: The magnetization histogram (squares and lines) for the 2D Ising model with J = 1, B = 0.01, and kT = 2.5 on a 
square lattice of size: (a) L — 10; (b) L — 20. The crosses show the prediction of the fluctuation relation (|115|) . 



In Fig. [TUl the histogram of the magnetization is shown at kT = 2.5 close to the transition. We observe that the 
distribution is still bimodal for the small system, but it becomes unimodal as the size increases. 
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FIG. 11: The magnetization histogram (squares and lines) for the 2D Ising model with J = 1, B = 0.01, and kT — 3 on a 
square lattice of size: (a) L — 10; (b) L — 20. The crosses show the prediction of the fluctuation relation (|115|) . 



Finally, Fig. [TT] shows the magnetization histogram at kT = 3 above the transition where the bimodality has 
completely disappeared. 

Away from the critical temperature, the peaks of the magnetization distribution become sharper as the size increases, 
which is expected from the large-deviation expression (|28p . 

In all cases, the magnetization distribution computed by Monte-Carlo algorithm is compared with the prediction 
of the fluctuation relation (|2"B1) according to which 

p B (m) = ps(-m) exp(2f3BNm) (115) 

The nice agreement observed in Figs. CTTTl confirms the validity of the fluctuation relation in the 2D Ising model. 



IX. CONCLUSIONS 



In this paper, relationships have been established that express the breaking of discrete symmetris by an external field 
in equilibrium statistical mechanics. To be specific, we have considered spin systems described by Gibbsian canonical 
equilibrium states based on a Hamiltonian function of the form given by Eq. (J2J) and satisfying the symmetry conditions 
(01- © under spin reversal 

With these assumptions, the exact fluctuation relation (|26|) is deduced for the probability distribution of the 
magnetization. In the large-system limit, this fluctuation relation implies symmetry relations for the large-deviation 
function and the cumulant generating function of the fluctuating magnetization per spin. 

Furthermore, a concept of coentropy or spin-reversed entropy is introduced, which forms a Kullback-Leibler diver- 
gence if combined with the standard thermodynamic entropy. Under the assumptions considered, this Kullback-Leibler 
divergence is related to the product of the external field with the average magnetization. By the non-negativity of the 
Kullback-Leibler divergence, this product is thus always non negative, implying a paramagnetic response to an exter- 
nal field. While the entropy is a measure of disorder in typical spin configurations, the coentropy characterizes the 
rarity of the corresponding spin-reversed configurations. The larger the coentropy, the rarer the spin configurations 
opposite to the most probable ones. 

These results have been applied to four different spin models in the presence of an external magnetic field breaking 
the spin-reversal symmetry: the non-interacting spin model, the ID and 2D Ising models, as well as the Curie- Weiss 
model. These last two models feature a paramagnetic-ferromagnetic phase transition. In the ferromagnetic phase, 
the magnetization distribution is bimodal, but the fluctuation relation is always satisfied. In this way, the study 
of these models show how the established relationships characterize the breaking of spin-reversal symmetry by the 
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external magnetic field. In spite of the symmetry breaking phenomenon, the relationships show that the underlying 
fundamental symmetry ((5]) continues to manifest itself. 





Equilibrium 


Nonequilibrium 


broken symmetry 
external control parameter 
order parameter 

non-negative quantity 

fluctuation relation 


spin reversal 
external magnetic field B 
average magnetization (m)s 

i( S R -s) = 2/3B(m) s 
P B {M) = P B {-M)e 2fiBM 


time reversal 
affinity A 
average current (J) a 

entropy production \ = h R — h = A(J) a 

P A (J) = PA(-J)e AJ 



TABLE I: Comparison between the equilibrium and nonequilibrium broken symmetries, k denotes Boltzmann's constant and 
ft = (fcT) -1 the inverse temperature. 



As shown in Table HI there is an analogy with similar results already known in nonequilibrium statistical mechanics 
In this other context, fluctuation relations have been obtained for currents flowing across open systems in 
nonequilibrium steady states breaking the time-reversal symmetry. In the analogy, the magnetization corresponds 
to the current and the external magnetic field to the thermodynamic force or affinity inducing non-zero values of 
the average current. In this context as well, a dynamical coentropy or time-reversed entropy per unit time has been 
previously introduced giving the thermodynamic entropy production as a Kullback-Leibler divergence Q. 

This analogy emphasizes the role of these relationships in characterizing the breaking of some discrete Z2 symmetry 
at the statistical level of description. In this regard, we may expect that such relationships can be used to study 
other discrete Z2 symmetries, in particular, at equilibrium. Instead of systems with magnetic moments, systems with 
electric dipoles featuring transitions towards a ferroelectric phase could also be considered. Moreover, the present 
relationships can be directly extended to lattice gases where spin reversal corresponds to the exchange of empty and 
occupied sites (l9| . Further discrete symmetries can also be envisaged such as molecular chirality in order to investigate 
the problem of the emergence of homochirality in physico-chemical systems, or matter-antimatter asymmetry under 
early cosmological conditions. The present results also suggest to study combined broken discrete symmetries such 
as spin and time reversals in nonequilibrium steady states with an external magnetic field. Beyond these contexts, 
similar symmetry relations have already been obtained for multifractal spectra at random critical points 23]. All 
these advances open new perspectives for a better understanding of symmetry breaking phenomena. 
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